Solitary Wave Interactions In Dispersive Equations Using Manton's Approach 
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We generalize the approach first proposed by Manton [Nuc. Phys. B 150, 397 (1979)] to 
compute solitary wave interactions in translationally invariant, dispersive equations that support 
such localized solutions. The approach is illustrated using as examples solitons in the Korteweg-de 
Vries equation, standing waves in the nonlinear Schrodinger equation and kinks as well as breathers 
of the sine-Gordon equation. 



Introduction. Dispersive wave nonlinear partial differ- 
ential equations (PDEs) describe a variety of physical 
systems in atomic, optical, molecular, solid state and 
wave physics as well as in fluid dynamics, biophysics, 
plasma physics, high energy physics and astrophysics 
among others [1,2]. Often in these settings, it is of par- 
ticular interest to examine the dynamics and interactions 
of spatially localized, and possibly travelling in time, so- 
lutions which can represent bits of information, moving 
Bose-Einstein condensates, elementary particles or water 
waves [3,4]. Particularly the interactions between the 
solitary waves are especially important, since, e.g., in 
optical communications avoiding such interactions may 
reduce the bit error rate [5]; in Bose-Einstein conden- 
sates such interactions change significantly the form of 
the wavefunction [6], while in high energy physics mod- 
els, the interactions are used to monitor the collisions of 
elementary particles [7] . 

There is a vast amount of literature regarding the in- 
teractions and collisions of the above mentioned solitary 
waves. Many reviews of continuum [8,9] and discrete [10] 
systems contain some of this enormous volume of work. 
The techniques that have been used to identify the na- 
ture of such interactions are also rather diverse ranging 
from perturbation theoretic ones as, e.g., in the work of 
Ref. [11], to variational ones [9,12,13], to more rigorous 
ones, using the Fredholm alternative [14] or Lin's method 
as in Ref. [15]. Typically, these interactions asymptoti- 
cally follow the tails of the waves, which in most cases are 
exponential. This, in turn, results in writing down Toda 
lattice type equations at the "mesoscopic" level for lat- 
tices of coherent nonlinear waves, see, e.g., Refs. [13,16] 
and references therein. 

While the solitary wave interactions have been exten- 
sively studied in the past, in the present communication 
we aim at presenting a different viewpoint on this topic. 
Our aim, in particular, is to explicitly focus on calcu- 
lating the tail-tail interactions between the waves, using 
the method proposed by Manton in Ref. [17], based on 
the earlier work of Refs. [18-20] and systematize it as a 
general method that can be straightforwardly applied to 
any nonlinear dispersive wave equation that has a num- 
ber of characteristics (which will be analyzed/explained). 



We believe that this method provides a very simple, yet 
elegant and useful tool that can be generically used in 
this large class of models and hence would be of value to 
researchers in a variety of disciplines. 

The structure of our presentation is as follows: initially 
we repeat Manton's formulation for the kink-antikink in- 
teraction in the sine-Gordon equation, highlighting some 
key and subtle points. We then proceed to an analogous 
calculation for the case of solitons in the Korteweg-de 
Vries (KdV) equation. We then move to the realm of 
breathers starting with the standing waves (with triv- 
ial periodicity) in the case of the nonlinear- Schrodinger 
(NLS) equation. Thereafter, we study the interaction of 
genuinely breathing structures such as the breathers of 
the sine-Gordon equation. All of our results are corrob- 
orated with numerical simulations. Finally, we conclude 
our presentation with a summary of main findings and 
some intriguing questions for future study. 

Manton's Formulation for sine-Gordon Kinks. At the 
heart of Manton's calculation of the interaction energy is 
the use of the definition of the linear momentum of the 
wave equation at hand. Computing the time derivative of 
the momentum P for an interval containing one solitary 
wave, and deducing the contribution to it from the sec- 
ond wave, we can infer the force exerted on the soliton 
from its neighbor. In particular, for Lorentz-invariant 
equations (e.g. Klein-Gordon equations) of the form: 



u t t 



V'(u), 



the linear momentum can be written as: 



P = — / u t u x dx. 



(1) 



(2) 



The total integral along the line is a conserved quantity 
(when the model is translationally invariant, an assump- 
tion necessary for this approach to work). Consider a 
soliton centered around £ = (u^'), and an antisoliton 
centered at £ = A£ (it*- 2 -*), then using an interval (a, b) 
such that aCO and <C b <C A£, we find that in that 
interval (as shown in Ref. [17]): 



dP 



(3) 
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If we then take into account the asymptotic form of the 
kink-like waves, which for V(u) — 1 — cos(u) are of the 
form u = 4arctan(e ±2: ), we find that uW « 4exp(— x) 
and ~ 4exp(— (x — A£)) in the region between the 
solitons (and hence at x — b). At x = a the contribution 
will be null as a — > — oo. Manton's end result is thus: 



UJ 1 



dP 

~dt 



32exp(-A£), 



(4) 



(which is independent of b) and hence the potential of 
interaction is V(A£) — — 32exp(— A£). 

There are some subtleties here that we would like to 
highlight: 

• Notice that dP/dt = dV/dAt; has been used in the 
derivation of the potential, as opposed to the usual (— ) 
sign on the right hand side. This is because, the way the 
momentum is defined, dP/dt > which means that the 
momentum at x — b is increasing, which in turn means 
that the solitary wave is approaching x — b, because of 
the interaction and hence A£ is decreasing (i.e., when 
dP/dt > 0, the resulting force in the dynamics of A£ is 
negative). This sign change accounts for the above for- 
mula used in the derivation of the potential. This subtle 
point should always be taken into consideration when us- 
ing Manton's method. 

• Secondly, if we are to examine the dynamics of A£, 
then the corresponding Newton equation should account 
not only for the potential, but also for the "mass" 
of the solitary wave. In the case at hand (kinks of 
the sine-Gordon equation), the mass M = J u 2 dx = 

I [hi + (1 " cos(u))] is 8. 

• Finally, for the dynamical equation of A£, we should 
take into consideration that we have only computed the 
force on the "left soliton" from the right one. However, 
there is an equal and opposite force on the right solitary 
wave and hence we should use a factor of 2 when writing 
the equation for the acceleration of A£. 
Incorporating all these points in the equation for the sep- 
aration between the solitary waves we obtain that, Man- 
ton's formulation yields in general: 



A£ = - 



2 dP 
M~dt' 



(5) 



where the overdot denotes temporal derivative. In the 
particular case of sine-Gordon kink-antikink interaction, 
Eq. (5) yields 



A£ = -8exp(-A£) 



(6) 



A numerical example illustrating/validating the useful- 
ness of this formula is given in Fig. 1. It is worth remark- 
ing here that for all Lorentz-invariant theories, while the 
kink-antikink interaction is attractive, the kink-kink in- 
teraction is always repulsive. 




FIG. 1. Kink-antikink attractive interaction in the 
sine-Gordon equation. The kinks are initialized at a distance 
of ~ 10.4 and are allowed to interact. The top panel shows 
their separation (numerically measured using a local center 
of mass approach weighted on the energy density) as a func- 
tion of time. The solid line shows the numerical result, while 
the dashed one indicates the result of integrating Eq. (6). 
The bottom panel shows the space-time (x-t) contour plot of 
the energy density, indicating the tendency of the waves to 
approach each other. 

Soliton Interaction in the KdV Equation. Another in- 
teresting example that illustrates some additional sub- 
tleties of the Manton approach concerns the solitary wave 
interactions in the KdV equation (most often associated 
with water waves [3,4]). Following the same path as 
above, the definition of momentum in this case is given 
by P — J u 2 dx, hence we compute dP/dt. Given the 
form of KdV (in travelling wave frame with velocity C): 



Ui 

we find that 
dP rb 
~dt ~ 



6uu x + Cu x 



uu t 



[Au 3 - 2uu x 



Cu l 



(7) 



(8) 



For a two-soliton decomposition u — u^> + and ne- 
glecting (exponentially smaller) higher order terms, we 
obtain the expression: 

rlP 

^- w 2«W«W - 2( U «4 2 j + u^ui 1 }) + 2Cu^u^ . (9) 

Finally, given the form of the soliton u = 
(C/2)sech 2 (VCx/2), the expression of Eq. (9) is evalu- 
ated as: 



dP / 

_ = -16C 3 exp (-VCA£ 



(10) 



Using the expression of Eq. (5), we would be immedi- 
ately inclined to write (given that the mass of the soliton 
is J udx — 2\fC) 

At = 16C 5/2 cxp f-VC A^j . (11) 

However, we argue that when following Manton's formal- 
ism, one should respect Ehrenfest's theorem. In particu- 
lar, in the case of KdV, it is true that: 



2 



dt 



idx — 3P 



dt 2 



xudx = 3' 



dP 
dt' 



(12) 



The left hand side indeed evaluates to MA£ here (as 
well as for sine-Gordon), however the right hand side 
should be multiplied by a factor of 3. Hence, the result- 
ing dynamical equation for the evolution of the soliton 
displacement should read: 



A£ = 48C 5/2 exp (-VCA^j 



(13) 



A numerical example of the dynamics of KdV solitons is 
given in Fig. 2. 



waves, and the relative phase (f> between them has been 
incorporated in U2- Then, one obtains: 



dP 

— =87 7 4 exp(- ?7 AO, 



(17) 



which results in the dynamical equation for the separa- 
tion (using Eq. (5)) of the form: 



A£ = -8?y 3 exp (—rjA£) cos < 



(18) 



This equation is identical to the one obtained in [21] 
through variational and perturbation techniques [cf. Eq. 
(5) of [21]]. 

We have also examined numerically the validity of the 
results. A typical computation is summarized in Fig. 3. 
Notice that in the case of solitons in phase, the interac- 
tion is attractive while it is repulsive when the solitons 
are (it) out of phase. This is also observed in the numer- 
ical experiments. 




FIG. 2. The (repulsive) two-soliton interaction in the KdV 
equation. The two panels are similar to those of Fig. 1. The 
small jumps in the top panel are due to the discretization 
used to emulate the continuum equation. The two solitons 
are initialized at a distance of A£ = 9. The bottom panel 
shows the space-time (x-t) contour plot of the field it. 



Standing Wave Interaction in the NLS Equation. In 
the case of the NLS equation (a model relevant predom- 
inantly to nonlinear optics and plasma physics [5], but 
also more recently to atomic physics as well [2] ) , the role 
of the mass is played by the squared L 2 norm, hence mass 
and momentum are defined as follows: 

M = J \u\ 2 dx; P = % - J {uu% - u*u x )dx, (14) 



for the dynamical equation: 



1 



lUt 



--u TT . - u u. 



(15) 



We can then compute similarly as above (up to higher 
order exponentially smaller terms): 

= \ [«<x + u xx u*-2\u x \ 2 ] b a . (16) 

We again use the standard 2-soliton decomposition, u — 
ir 1 ) + where = i]sec}i(rix)exp(ir) 2 t/2) 1 = 
■qsech(ri(x — A£)) exp(in 2 t/2) exp(i(j)) are the standing 




FIG. 3. Same as in Fig. 1, but for the (in phase) stand- 
ing waves of the NLS equation. The latter were initialized 
at a distance of 6. The bottom panel shows the space-time 
contour plots of the square modulus of the wavefunction u. 



Breather Interaction in the sine-Gordon Equation. In 
the above examples we established cases where the inter- 
action between the solitary waves was previously known 
and used them to analyze some of the key technical points 
of the Manton approach. We now turn to an example 
that, to the best of our knowledge, has not been treated 
analytically before and which concerns, in particular, the 
breather-breather interaction in the sine-Gordon equa- 
tion. Such solutions of Eq. (1) are given by the form: 



u = 4 arctan 



j£)sech y\f 



1 



(19) 



where u> is the breather frequency. A fundamental differ- 
ence in this case is that the solutions are genuinely time 
dependent. However, the Manton formulation can still 
be carried through with the corresponding momentum 
derivative evaluated as: 



dP 
~dt 



,(!)„, ( 2 ) 



,(1)„(2) 



+ {u 



(1) 

XX 



)u 2 



(20) 
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result of the decomposition 



As 

+ and the asymptotics vP^ 

Sy/1 — uj 2 sin(wt) exp (— \/l — lu 2 x) /uj and w 
8Vl — w 2 sin(wi) exp (— yl — lu 2 (x — A£)) /uj, we ob- 
tain: 



dP 64(1 -w 2 ) 5 



1 



I - --^ cos(2wi) ) 
1 — u> z ) 

exp(-v / l - uj 2 AQ. 



(21) 



From this calculation, we can infer that while the poten- 
tial of interaction between two sine-Gordon breathers is 
non-autonomous, it has a well defined attractive (when 
the breathers are in phase) average of V — — (64/w)(l — 
uj 2 ) 3 ' 2 exp(— Vl — w 2 A£). Furthermore, using Eq. (5) 
and the breather mass M = 16^ — uj 2 , we can infer the 
dynamical equation governing the inter-breather separa- 
tion as: 



8(1 -uj 2 ) 3 / 2 ( 1 .„ , N 

■ ^ (l- T — rCos(2^ 

x exp (-- v/l - w 2 A£). (22) 



A numerical experiment illustrating the comparison of 
Eq. (22) with the PDE result is shown in Fig. 4. 




FIG. 4. Same as in Fig. 1, but for the (attractive) in- 
teraction of in-phase breathers of the sine-Gordon equation. 
The breathers are initialized at a distance of ~ 16.45. The 
dashed line is the theoretical prediction of Eq. (22). The small 
jumps in the numerical result are an artifact of the discretiza- 
tion used to approximate the continuum equation. However, 
it can be clearly seen that the theoretical prediction follows 
well the full numerical evolution at least for times t < 35, 
where the breathers are sufficiently well separated that the 
interaction has not changed their profile significantly (and for 
which they maintain their individual character). 



Conclusions. In this short communication, we revisited 
the topic of solitary wave interactions for exponentially 
localized solutions (even though that is not absolutely 
necessary; however, it is typically the case) of transla- 
tionally invariant, nonlinear dispersive wave equations. 
We showed how to systematically exploit the momentum 



integral to find the force exerted on one of the waves by 
the other and how to establish the dynamical equation 
of motion of the inter-soliton distance, using the mass 
integral and Ehrenfest's theorem [notice that due to the 
structure of the energy-momentum tensor for Lorentz in- 
variant equations, the analog of such a theorem in the 
latter case is given by d/dt J x£dx = P, where £ is the 
energy density]. We demonstrated the use of the ap- 
proach in a number of well-established situations, includ- 
ing the kink interaction in the sine-Gordon equation, the 
soliton interaction in the Korteweg-de Vries equation and 
the standing wave dynamics in the nonlinear Schrodinger 
equation. The method was also applied to obtain an an- 
alytical expression for the breather-breather interaction 
in the sine-Gordon model. 

Naturally, the approach has a number of limitations: 
for instance it cannot be directly applied in the pres- 
ence of spatially dependent potentials. Similarly, it can- 
not be straightforwardly implemented in cases where 
Ehrenfest's theorem (or Lorentz invariance) cannot be 
established. A prominent such example is, for in- 
stance, the modified Korteweg-de Vries equation wherein 
d/dt J xudx = 2 J u 3 dx, whereas the momentum is given 
by J u 2 dx. Certainly, it is also of interest to generalize 
the approach to multiple dimensions in the spirit of Ref. 
[12]. Finally, an alternative method to compute the inter- 
action energy of solitary waves is through direct energy 
arguments: in particular, if the corresponding soliton lat- 
tice solution exists, one can compute the energy of such 
a solution, inside a period. The leading term will then be 
the energy of a single soliton, while the corrections will 
correspond to the soliton interaction energy. It would 
be interesting to compare/validate the results of such a 
method with those of Manton's formalism. These topics 
will be deferred to future publications. 
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